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ABSTRACT 


Detection of edges under noisy environments has been gaining lot of 
prominence in the recent past in most of the image and video processing 
applications. In this work a novel approach based on the distribution of 
intensity values and their corresponding positions has been proposed for 
distinguishing the edge pixels from the grey scale images. Separate 


histogram has been maintained for X and Y coordinates. The first order 

derivative is applied over these histograms to distinguish the edge pixels. The 
Keywords: pixel with gradient distribution below a specific threshold value is selected as 
an edge pixel. This method is found to work well in case of both noiseless 
and noisy images. Hence this method is able to perceive the underlying 
information in case of noisy images. The proposed algorithm can be used for 
both low and high resolution images. However, the performance of the 
algorithm is more evident in high resolution image. A general analysis of the 
proposed method has been conducted for arbitrary images. The major 
application of the proposed work can be used for the applications that doesn’t 
need any preprocessing or to avoid any loss of information like in medical 
image analysis as it contemplate towards every intensity bin to trace the 
edges present in the histogram of the image rather than the overall image 
concerning for direct edge tracing. The results have been compared with 
canny algorithm which is most commonly used for edge detection. 
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1, INTRODUCTION 

Edges play a dominant role in the image processing and computer vision field. It is used in the 
object detection, shape detection and more. There are algorithms still in use for better edge detection. Edges 
within the noisy image are marked along with non-salient information. Image acquisition at some typical 
situation leads to noises. The image may be captured in the smoke or fog environment. The semantics part of 
the image may lose during its edge detection. In this paper, a new algorithm is developed based on the 
histogram technique. The result is based on the presence of intensity distribution, even the image is hazy. 
Even human perception may fail at just looking upon the image. The boundary of the object may be difficult 
to perceive due to heavy fog or smog. The computational approach of this algorithm determines the edges in 
the hazy situation and gives good results. There are results which have been discussed here. The astronomical 
image edge detection gives a great result compared with existing model. Texture determination is also 
difficult in noisy images. Even though there are few methods which are implemented to extract the edges, the 
proposed approach besides gives a better result. The texture content is preserved during this edge extraction 
method. At certain situation where the texture content is more important, the edge will be determined even in 
the obscure condition. Edge information is lost in low-resolution images when the images are converted from 
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higher-resolution. The proposed approach gives a far better result for the edge enriched images than the 
existing methodology. 


1.1 Background 

Guo, et.al [1], proposed the method to detect shadows and its removal from the still images. 
Illumination discontinuities are not considered for the shadow region. The comparison is performed if the 
segmented region and the others are having similar pattern region. Soft matting is applied here. 

Gupta and Das [2] proposed a textured edge detection algorithm using a bank of multi resolution 
filters and self-organizing map (SOM). 2D filtering traces the possible edges in all directions. SOM performs 
the dimensionality reduction by producing a 1-D map which groups the similar items together. Gabor filter 
and DWT is used as an advantage for texture edge detection. This produces accurate results of texture edges 
for natural texture images. 

Huanhuan and Nenghai [3] presented a fast algorithm to find the edge saliency maps for a natural 
image. This approach integrates edge gradient value, edge segment length and edge density. They claim that 
the approach outperforms other edge saliency detection methods. This gives more complete object edges and 
also gives the non-salient edges. 

Krishnamurthy, et.al [4] presented a new edge detection algorithm for oceanographic features 
present in infrared (IR) images obtained from the Advanced Very High Resolution Radiometer (AVHRR). 
The Histogram Based Morphological Edge Detector (HMED), extracts all the weak gradients which are not 
distinguished by the conventional edge detectors. This extracts the edge fine structures without introducing 
any spurious edge pixels. 

Kim and Park [5] presented an efficient video segmentation strategy based on new edge features. 
The introduced edge information is determined from the analysis of the local histogram. According to the 
authors, the edge feature is extracted where a change of intensity or color is weak. There exists an advantage 
of segmenting fewer regions compared with the other existing techniques. 

Xu, et.al [6] presented a unique way of detecting the edges. Here the canny edge detection algorithm 
is applied to every individual block instead for the entire image. The adaptive threshold selection method 
determines the high and low thresholds while processing for every individual blocks. A set of work [14], [12] 
and [13] has been reported based on Deriche filters that have been derived using Canny’s criteria and 
implemented on ASIC-based platforms. R. Deriche [14] presented a recursive algorithm for edge detection 
and from Canny’s design [18] a solution for an infinite extent filter leads to an optimal operator in one 
dimension and a two recursive filters 1s used for implementation is shown. 

L. Torres, et.al [12] presented a design of a real time image processing circuit based on an optimized 
Canny-Deriche filter for ramp edge detection. Retiming method is used in [12] for achieving very high speed 
filtering. Deriche uses the recursive filter to reduce the computation cost. F.G.Lorca, et.al [13] optimised the 
Deriche filter both algorithmic and architectural aspects and organises the filter at the 2D and 1D levels. 

Carrato and Ramponi [14] proposed a novel approach to recognise the edges based on the evaluation 
of the proportion of the data present on a mask which scans the entire image. The generalised higher order 
statistical moment is used to search the zero crossings for presenting the edges of the original image. 

Jain, et.al [12] presented a novel idea to detect the foreground objects using sub pixel edge maps on 
a video sequence. Sub pixel edge models are traced efficiently using a modified canny edge detector as it 
gives a better result in case of sudden illumination changes. 

Lee, et.al [13] proposes a new unsupervised segmentation method for hyper spectral images based 
on dominant edge information. The principal component analysis is applied first and the dominant Eigen 
images are selected. The edge operators and histogram equalizer are applied to the selected Eigen images, 
which produce edge images. These edge images are combined, a binary image is obtained. Morphological 
Operators are applied to these binary images to remove erroneous edges. With the edge information of the 
dominant Eigen images, segmentation is made. 

Leone and Distante [14], presents a new approach for automatic segmentation of potential shadow 
points in the visual surveillance environment by taking the background difference into considerations, for all 
moving pixels. Shadows are described here as half transparent region of the background textural 
characteristics and also claims that the algorithm is unaffected by scene type or any illuminations. 

Vidya, et.al [15] emphasize the edge detection using masks on hexagonal grid (using spiral 
addressing scheme) and gabor wavelet based edge detection is proposed for the hexagonal grid. Abhishek, 
et.al [16] presented a new approach for extracting the edge information. This paper has a detailed study about 
the proposed intelligent algorithm. Sowmya, et.al [17] presented an application of sparse banded filter 
matrices in edge detection. The filter design is formulated in terms of banded matrices. The continuous edges 
are attained without any parameter tuning. 
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The most popular technique used for edge detection is canny edge detection method. This gives low 
error rate as proposed by Canny [18]. In this method, Image is smoothened by a Gaussian filter and the 
intensity gradient along x direction and y direction is computed using Sobel kernel. Gradient is the change in 
brightness over a series of pixels. Direction can be found using the gradient along x and y directions. Non 
maxima suppression is performed by removing pixels that are not considered to be part of an edge. A binary 
image is formed with the thin edge. The thin edge is maintained by removing the edge pixels that are not 
along the edge direction. Canny preserves the dominant edges. Canny detection methodology uses the 
gradient values in determining edges, where some low level gradients will be failed to recognize. Due to the 
loss of quality during the acquisition of the image, distortions of the intensity will be more that will make the 
edge detection worse. 

The major advantage of the proposed work is the presence of intensity distribution over the image. 
The low level gradients are identified easily as the algorithm performs computation for all intensities 
individually using histogram technique. 


1.2 Problem 

The focal length of a camera affects the image in case of locating the far away objects. Edge 
detection over such images gives poor results. The images affected by noises will hold its most anticipated 
region less obscure. The proposed method recovers edge information in these noisy images. This method is 
found to be less affected by the illumination. 


1.3. Proposed Solution 

The proposed method can define the semantics part clearly along with the edges. Histogram 
technique is used to give the importance to all intensities. Hence, the weak gradients are also identified along 
with the salient edges. The nearer coherent intensity blocks sometimes deceived as separate regions due to 
the presence of noises or poor quality of the image. The proposed approach founds the nearer coherent region 
to avoid false noises and selects the better edge results. This is one of the major advantages over the proposed 
approach. This illusion of different region is clearly interpreted in Figure 8. These can be avoided by 
choosing better threshold value to trace out true edges. If threshold value is low then algorithm traces every 
nearer coherent intensity blocks boundary as edges and prediction of false noises will be more. Thus the 
threshold value is chosen bigger in order to avoid more false edge in the image. The canny edge detection 
algorithm works good but for certain noisy images, it gives the edge information for noises also. The 
performance of this approach is more evident when considering high resolution images. 

The edge in the foreground shows much dominance than the background. Thus edges of the 
foreground objects will be traced evidently irrespective of the noises present in the image. Diagonal edges 
can also be traced using this proposed algorithm and provides obscure edge information as well. The dark 
textured background affected by shadows with and without noises can be traced. Experimental results shows 
that for lower resolution images, the parameters used for comparison here gives lower values indicates false 
edges are more. But for high resolution images, result shows that values are nearer to canny and sobel 
approach for common edge pixel present in both vertical and horizontal edge information. So, these results 
evidently shows, traced edges are nearly closer to true edges. 


2. RESEARCH METHOD 

The proposed approach is tested with BSDS300 Berkeley segmentation dataset and also with self 
acquired dataset which was made for this research. The images with various backgrounds and its noise 
model, texture model, astronomical images are tested with the proposed algorithm. Some of the sample 
images used [25], [20] and [21] here are shown in Figure 1. 
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Figure |. (a) Alarm Clock, (b) 208001, (c) 160068, (d) 58060, (e) Cup, (f) Amazon Parrot, (g) Gorilla, (h) 
NGC 6543, (4) Orangutan, (j) Monkey, (k) NGC 6384-Spiral Galaxy. 


The gray scale image is taken as input by the algorithm. The edge pixels are selected based on the 
localization of the intensities. The proposed technique will deliver two types of output: vertical edges and 
horizontal edges. Combination of these edges results in new kind of edge information. Results are unique. 
The overall system design is shown in Figure 2. 





Figure 2. System Design. 


The proposed method uses histogram for analyzing the intensity distribution. An edge has been 
determined from the intensity bin of the histogram rather generating directly from the image. Pixel 
coordinates are collected, for every intensity ranges. The distribution of these pixel values will not be 
continuous. This discontinuity determines the beginning and ending of edges all over the image in the 
proposed approach. Threshold value is finalized to avoid setting unnecessary pixels as edges. The image 1s 
recreated with edge distinguishing pixels from the histogram. 


Let S be the set of pixel locations for the intensity p, 
Sp = {X,, Xo, X3, aia Xn} 


i€1,2..m 
Where xX, = {Xi, yi} for \ €1,2,..07 
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Let N be the new image, 


T be the Threshold value 
_ §1,1 XR) — Xe) < T 
Ley) 0 , Otherwise 


The image N(x, y) gives the edge information for the given image. The threshold value T is selected 
to suppress the few local maximum of the gradient values, for every intensity range. /(X;,) and I(X;4, ) both 
should lie on the same row or same column. 

Steps involved in edge distinguishing algorithm are as follows: 
1. Check the intensity values over the entire image. 
2. If the values stored in histogram are consecutive and distributed with minimum gradient, set pixel 
value as black otherwise change to white. 
3. Iterate through step 2 until last intensity value of the histogram bin. 

The threshold value is chosen within the range of the first order derivative. It can be selected 
either as a local threshold for the individual intensity or as a global threshold for all the intensity. Threshold 
is performed to avoid certain local maximum of the first order derivative. The first order derivative of the cup 
image is shown in Figure 4(g). The results can be improved by adjusting the threshold value. The flow chart 
for the algorithm is shown in Figure 3. There are some studies made on the directional derivative as in [22] 
that also gives comparatively good results. The fractional differential gradient operator is also used for 


extracting the edges [25]. 
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Figure 3. Flow chart of the proposed algorithm 
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Figure 4. (a) Original Image Cup (3264x1836), (b) Figure 5. (a) Original Image Clock (1544x868), (b) 


Horizontal Edges, (c) X coordinate for intensity Horizontal Edges, (c) X coordinate for intensity 
value 175 of 3 (a), (d) Zoomed view of 4 (c), (e) Y value 149 of 5 (a), (d) Zoomed view of 5 (c), (e) Y 
coordinate for intensity value 175 of 4 (a), (f) coordinate for intensity value 149 of 5 (a), (f) 
Zoomed view of 4 (e), (g) First order derivative of 4 Zoomed view of 5 (e), (g) First order derivative of 5 
(e), (h) zoomed view of 4 (g). (e), (h) zoomed view of 5 (g). 


The rejection of pixels will follow the conditions as follows: 

The first order derivative of the clock image is shown in Figure 5(g). The threshold value has been 
taken as 500 by brute force method and this won’t reject the pixels that has first order derivative less than the 
threshold value. The pixel intensity localization is considered if it is distributed continuously along the same 
row or column then the pixel is considered for comparing the first order derivative otherwise the pixel is 
rejected as non-salient edge. 

Consider the sample intensity values of an Image as in Figure 6. The proposed algorithm is 
implemented in these values. The algorithm checks the consecutiveness of these values. If the intensity 
values are stored in consecutive locations, set pixel value as white at initial and terminal else set as black. 
Consider the intensity value 200, set the initial and terminal as 255 and remaining as 0. This sample value is 
verified, the result of the processed image is shown in Figure 7. The edge pixels are marked in white color. 






cS 


}255 |o | 25s | 255 | 255 |o | 255 
j25s Jo 25s [255 | 255 |o [25s 











Figure 6. Sample values and its processed output. Figure 7. Sample Image and its processed 
output. 
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The alarm clock image is taken for interpreting the pixel selection. The intensity value 149 is 
highlighted in Figure 8 (a). Figure 8 (a) and Figure 8 (b) show the comparison about the disappearances of 
certain pixels. Most of the pixels are not distributed continuously and they are rejected. The pixel that 
satisfies the constraint is marked in the resultant image is shown in Figure 8 (d). 

This algorithm is tested with self-acquired data set. The gray level value of the alarm clock is given 
as input. It gives the result as shown in Figure 10. This gives vertical edges and horizontal edges. This edge 
information is utilized for developing new results by combining both the edges. It gives collective results as 
shown in Figure 10 (d). Comparison is made with the existent methodology is shown in Figure 10. Results 
are capable of recovering the trace of shadow, fine than the existing methodology. 





8(b) 





8(c) 8(d) 9(c) 


Figure 8. (a) Intensity value 149 of alarm clock Figure 9. (a) Zoomed view of the alarm clock image 
image, (b) Selected edge distinguished pixels, (c) (b) vertical edges for 9 (a) zoomed region (c) 
Region selected for zooming from 8(b), (d) Zoomed Horizontal edge for 9 (a) zoomed region (d) 
view of the selected region in 8(c). Combined results of both the vertical and horizontal 
edges. 





10(e) 10(f) 


Figure 10. (a) Original Image Clock (1544x868), (b) Vertical Edges, (c) Horizontal Edges, (d) Combination 
of Vertical Edges and Horizontal Edges, (e) Canny Edges, (f) Sobel Edges. 


3. RESULTS AND ANALYSIS 

The algorithm is tested with standard dataset and also with the self-acquired dataset. The high 
resolution images are giving a better result. The results are good for higher resolution images since there is 
no loss of information or the original image is not under sampled. The intensity values are directly processed 
for finding its distribution to determine the edges. The horizontal and vertical edges are providing unique 
information. Few images gave a better detail in vertical edge information. Horizontal edge information had 


Edge Detection Using Histogram Localization (Mohamed Asharudeen J) 


348 O ISSN: 2502-4752 


given enhanced results for other images. In this research the tests are done with different types of images. 
The result of the proposed algorithm is also utilized to bring two hybrid results. Combining both horizontal 
and vertical edge information had given stronger edge details. Common pixels present in both vertical edge 
information and horizontal edge information gives weaker edge details. For certain complicated images, the 
common pixel edge information results are clear compared to the other results. 

Shadows are dark transparent part of the background present within the image. Shadow tracing of an 
object as shown in Figure 10 and Figure 11 is also possible with the algorithm. When compared the results 
with the existing methodology, shadow edge tracing gave fine results in this algorithm. The images shown in 
Figure 10 and Figure 11 are rich in illuminations. The edge detection for these images with the existing 
approach traces the unwanted edges along with salient information. Even though the canny approach 
performs Gaussian blurring [24] over the input image for edge detection, the edge traced by the proposed 
algorithm is improved for higher resolution images without implementing any preprocessing technique over 
the image. 

Gorilla image as shown in Figure 17 had many noises around its body. Due to the intensity 
distribution of the hair, the output showed the intensity variations of every hair. The proposed approach 
brings the detail of the depth information in the image which is not shown in the existing model. This detail 
of the depth comparison is shown in Figure 17 with red markings. The comparison over the existing model is 
also shown for orangutan image in Figure 12 respectively. 





L1(a) Li(b} Li(c) 
9 





1l(e} Liff 





Figure 11. (a) Original Image Cup (3264x1836), (b) Vertical Edges, (c) Horizontal Edges, (d) Combination 
of Vertical Edges and Horizontal Edges, (e) Common pixels of Vertical Edges and Horizontal Edges, (f) 
Canny Edges, (g) Sobel Edges. 


The proposed algorithm is also tested with noisy images. The prominent information of the image is 
retained. The noisy images used for testing contained Gaussian noises and salt & pepper noises. The output 
consists of more noises due to the presence of larger variations in intensity as noise pixels. As stated earlier, 
for certain images as shown in Table 1, the combined results of vertical and horizontal edge information gave 
good result compared with the individual results. 

Texture contains more edge information than the typical region in an image. Texture tracing in this 
proposed approach gives good end result as traced in Table 2. The intensity distribution plays a dominant role 
in this method, even though the image is affected with noises. 


Indonesian J Elec Eng & Comp Sci, Vol. 11, No. 1, July 2018 : xx — xx 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 349 





i2{g) 


Figure 12. (a) Original Image Orangutan Figure 13. (a) Original Image NGC 6384-Spiral 
(1600x1200), (b) Vertical Edges, (c) Horizontal Galaxy (3871x1836), (b) Vertical Edges, (c) 
Edges, (d) Combination of Vertical Edges and Horizontal Edges, (d) Combination of Vertical Edges 
Horizontal Edges, (e) Common pixels of Vertical and Horizontal Edges, (e) Common pixels of 
Edges and Horizontal Edges, (f) Canny Edges, (g) Vertical Edges and Horizontal Edges, (f) Canny 
Sobel Edges. Edges, (g) Sobel Edges. 


Table 1. (a) Original Image (b) Vertical Edges, (c) Horizontal Edges, (d) Combination of Vertical Edges and 
Horizontal Edges, (e) Canny Edges, (f) Sobel Edges. 


Alarm Clock with Salt & Alarm Clock with Gaussian 


mavasel Pepper noise (1544x868) noise (1544x868) 


(a) 


(b) 


(C) 


(d) 


(e) 
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Table 2. Comparison of leopard image with noises for the existing and the proposed methodology, (a) 
Original Image, (b) Vertical Edges, (c) Sobel Edges, (d) Canny Edges. 
Dataset (a) (b) (c) (d) 
Leopard (481x321) | : _ 


Leopard with Gaussian noise 
(481x321) 


Leopard 
with Salt & Pepper noise 
(481x321) 





The proposed algorithm works well in the astronomical images also. The distant object captured 
using the space telescope is processed very well using image processing techniques. The complex processing 
of the astronomical data shows much space debris. The image shown in Figure 13 contains stars in the core 
part, and it is not easily detected with human perception sometimes. 





Late} Lait) L4tgi 


Figure 14. (a) Original Image NGC 6543 (720x780), Figure 15. (a) Cosmic dust hide stars near swirl of 


(b) Vertical Edges, (c) Horizontal Edges, (d) the NGC 6384-Spiral Galaxy shows with red 
Combination of Vertical Edges and Horizontal markings, (b) Combination of Vertical Edges and 
Edges, (e) Common pixels of Vertical Edges and Horizontal Edges shows hidden stars near swirl of 

Horizontal Edges, (f) Canny Edges, (g) Sobel Edges. the NGC 6384-Spiral Galaxy shows with red 


markings, (c) Canny Edges with partial results. 


The edge detection in these images brings a good vision of exploring the space objects. The core 
region of the spiral galaxy consists of stars that are not visible and it can be seen clearly using the proposed 
approach. Even the swirl in the core is not visible with direct human vision also it is not observable by the 
existing edge detection approach but the proposed algorithm shows clearly about swirls within the core of the 
image as shown in Figure 13 (d). If the astronomical data is in original resolution without any compression 
then we can see the hidden objects even better than the human perception. Some of the detected cosmic dust 
and stars are highlighted in red mark with comparison over the existing methodology and it is shown in 
Figure 15 and Figure 16 respectively. 
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Red marnings denote the wsibility of cosmic dust information with COMmparison Over The exehing 
and proposed marthoacds 


Red markings denote the visibility of neck and shoulder information with comparison over the 


existing and proposed methods 





Lola) 16{b) LGic] 17{a) 17 ic) 


Figure 16. (a) Original Image with obscure cosmic __ Figure 17. (a) Original Image with red marking about 
dust information in red marking, (b) Common pixels _ the shoulder and neck information, (b) Combination 


of Vertical Edges and Horizontal Edges with of Vertical Edges and Horizontal Edges with red 
obscure cosmic dust information is traced in red marking about the shoulder and neck information, (c) 
marking, (c) Canny Edge detection without obscure Canny Edges without the shoulder and neck 
cosmic dust information in red marking. information. 


Comparing the Lena image with different resolution in Table 3 with proposed and the existing 
approach, the higher resolution image has more prominent edge information. From top to bottom the image is 
filled with more noises given by the proposed approach. At the same time, when it is observed from bottom 
to top of the existing algorithm results, the results are better in low resolution images. 

The major intensity variation is clearly visible in the higher resolution image than the lower 
resolution image. The results of the existing method are noble at low resolution. The proposed results shows 
fine results at the high resolution images. There is much variation in intensity at high resolution images so the 
existing approaches are showing maximum number of unwanted edges. 


Table 3. Comparison of Lena image of different resolution with proposed and existing methodology, (a) 
Vertical Edges, (b) Horizontal Edges, (c) Combination of Vertical Edges and Horizontal Edges, (d) Common 
pixels of Vertical Edges and Horizontal Edges, (e) Canny Edges, (f) Sobel Edges. 


Dataset @) Fe © o 
Lena (256x256) Ff igaeeek L Pees Sh | aes. : 





Lena 
(1920x1920) 


Comparison of the performance of existing and the proposed edge detection methodology with 
respective ground truths are made and are shown in Table 4 and some of the image outputs are shown in 
Table 5. Ground truth images had been generated by using ImageJ tool. Berkeley Segmentation Dataset and 
few of its ground truth is taken for comparison as utilized in [25] by P. A. Khaire, et.al. Performance 
parameters used here are PSNR, Performance Ratio (PR) (Ratio of true to false edges), JACCARD and 
DICE. Performance results shown in Table 4 are taken for the low resolution images and high resolution 
images. PR is calculated as given in equation (1). True Edges (TE) referred here is edge pixels. False Edges 
(FE) referred here is the sum of non edge pixels identified as edges and edge pixels identified as non-edge 
pixels. JACCARD and DICE are calculated as given in equation (2) and (3). 
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TE 
PR =— xX 100 (1) 
FE 
Edge image AND ground truth image 
JACCARD = —————arrerowavrtrwr erat (2) 
Edge image OR ground truth image 
2(Edge image AND ground truth image 
DICE = cL Lea (3) 


Edge image + ground truth image 


The performance results of the proposed and the existing algorithm is visualized by projecting them 
in charts. There are four charts shown in Figure 18, Figure 19, Figure 20 and Figure 21 each shows how the 
proposed work is deviated from the existing approach. These results are acquired for test images of different 
resolutions. This also shows how the results are better from the existing approach and every result is not 
better in the proposed approach. The Figure 18 shows the PSNR value, if the error in the edge image is 
minimum, the PSNR value raises. The higher the PSNR values better the results. If the error is maximum, 
then the value becomes low. From Figure 18, it is obvious that the low resolution image gives low values as a 
result of more noises. The combination of both horizontal and vertical edges of lower resolution image gives 
negative values proves the presence of more noise. Mostly, common edge image of the proposed work has 
the plot line closer to the existing approach than the other results. Especially, higher resolution image have 
nearly closer to the existing plot line. The Figure 19 shows the PR value. 

The plot has some similar property as conferred in PSNR plot. The common edge image of higher 
resolution has the value closer to the existing approach. The plot in Figure 19 raises higher for certain images 
is actually depicting more false edges. The peak deviation is for the Lena image of dimension 256 X 256, as 
the proposed technique works better for high resolution images. The Figure 20 shows the JACCARD value. 
The Figure 21 shows the DICE value. The JACCARD index and DICE coefficient are similarity index 
describes how much the predicted image is similar to the ground truth image and its value lies between O and 
1. 





Figure 18. PSNR values for each test images. Figure 19. PR values for each test images. 


The common edge image line plot of higher resolution image mostly lies in between sobel and 
canny approach or nearer to the existing approach. The similar peak deviations in both of these plots are due 
to the result of the lower resolution lena image. From the overall performance, the common edge pixel of 
vertical and horizontal edges gives the most prominent edge information. For some cases the vertical edge or 
the horizontal edge gives better result depends on the applications. 





Figure 20. JACCARD values for each test images. Figure 21. DICE values for each test images. 
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The sample image taken for comparing the proposed method with canny edge detection result is 
shown in Figure 22 and Figure 23. The image has six parallel lines and another image is having an effect of 
salt and pepper noise. The major advantage over the proposed technique is that it rejects the false edges arise 
due to poor quality or the noises in the image. This can be explained by showing an intensity profile over the 
proposed output and the canny output over the test image. The spikes in the plot show how the edges are 
detected along with noises. Instead of six spikes for six parallel lines, the edge information for the noise 
image proves better than canny. The proposed approach is also tested with Gaussian noises in a similar way, 
but the results are not better as canny. 


22/(a) 





22(c} 23{c) 


23(d) 23(e) 
22d) 22[e) 


Figure 22. (a) Original Image (222 x 161) for Figure 23. (a) Original Image with salt & pepper 
comparing with canny, (b) Edge information of the noise (222 x 161) for comparing with canny, (b) 
proposed approach, (c) Edge information of the Edge information of the proposed approach, (c) Edge 
canny approach, (d) Intensity profile of the proposed information of the canny approach, (d) Intensity 
algorithm edge along 89" row, (e) Intensity profile of profile of the proposed algorithm edge along 89" 
the canny edge along 89" row. row, (e) Intensity profile of the canny edge along 89" 
row. 


Table 5. Some of the images used for comparing the proposed results with ground truth images along with 
the existing technique (a) Original Image, (b) Ground truth image, (c) Canny Edges, (d) Sobel Edges, (e) 
Vertical Edges, (f) Horizontal Edges, (g) Combination of Vertical Edges and Horizontal Edges, (h) Common 
pixels of Vertical Edges and Horizontal Edges. 


(C) (d) 





(1920 x 1920) 
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Table 4. Comparison of Performance Metrics with Existing Methodology 


Dataset Edge detection algorithm PSNR (in db) PR JACCARD DICE 
Church Canny 19.9279 15.2432 0.8566 0.9973 
Image (481 x Sobel 19.9181 11.1709 0.8995 0.9974 
321) Proposed method(Vertical edges) 4.6352 33.3313 0.7501 0.9978 
Proposed method(Horizontal edges) 6.8915 45.3854 0.6880 0.9970 
Proposed method(Vertical and Horizontal edges) 0.62953 34.9992 0.7407 0.9978 
Proposed method(Common edges found between 8.3499 44.1782 0.6940 0.9976 
Horizontal and Vertical edges) 

An Old man Canny 20.9899 11.5043 0.8768 0.9980 

Image (321 x Sobel 20.9797 8.9997 0.9174 0.9975 

481) Proposed method(Vertical edges) 4.1592 22.9026 0.8137 0.9984 
Proposed method(Horizontal edges) 3.2599 20.35 0.8307 0.9986 
Proposed method(Vertical and Horizontal edges) -1.2828 20.1888 0.8320 0.9987 
Proposed method(Common edges found between 6.3229 24.3949 0.8038 0.9981 
Horizontal and Vertical edges) 

Lena Image Canny 7.9271 6.5189 0.9035 0.9975 

(256 x 256) Sobel 7.9254 4.2425 0.9594 0.9988 
Proposed method(Vertical edges) 4.39 161.9216 0.3872 0.9929 
Proposed method(Horizontal edges) 3.4715 276.7512 0.2699 0.9927 
Proposed method(Vertical and Horizontal edges) -2.0542 700.2206 0.1310 0.9923 
Proposed method(Common edges found between 5.2869 91.3074 0.5267 0.9936 
Horizontal and Vertical edges) 

Lena Image Canny 18.1085 1553 0.9112 0.9955 

(1920 x Sobel 18.1086 1.5076 0.9865 0.9987 

1920) Proposed method(Vertical edges) 7.9197 24.7329 0.8028 0.9939 
Proposed method(Horizontal edges) 6.5812 35.4376 0.7395 0.9935 
Proposed method(Vertical and Horizontal edges) 2.7514 60.4318 0.6252 0.9929 
Proposed method(Common edges found between = 11.3259 9.0689 0.9174 0.9958 
Horizontal and Vertical edges) 

NGC 6543- Canny 14.5259 2.4963 0.9293 0.9963 

Spiral Galaxy Sobel 14.5253 2.2084 0.9785 0.9985 

Image Proposed method(Vertical edges) 8.9116 44.2423 0.7034 0.9932 

(720 x 780) Proposed method(Horizontal edges) 9.406 36.0305 0.7411 0.9935 
Proposed method(Vertical and Horizontal edges) 4.3749 73.9537 0.5893 0.9927 
Proposed method(Common edges found between 11.6363 16.977 0.8570 0.9947 
Horizontal and Vertical edges) 

NGC 6384- Canny 15.093 0.68035 0.8537 0.9953 

Spiral Galaxy Sobel 15.1 2.0564 0.9800 0.9989 

Image (3871 Proposed method(Vertical edges) 8.267 Da222 0.8115 0.9949 

x 1836) Proposed method(Horizontal edges) 9.2438 17.0705 0.8540 0.9955 
Proposed method(Vertical and Horizontal edges) 4.5697 38.4121 0.7233 0.9941 
Proposed method(Common edges found between = 12.529 6.0038 0.9424 0.9975 
Horizontal and Vertical edges) 

Clock Image Canny 20.5911 0.58404 0.9353 0.9954 

(1544 x 868) — Sobel 20.5969 1.0432 0.9897 0.9986 
Proposed method(Vertical edges) 9.017 17.4658 0.8480 0.9937 
Proposed method(Horizontal edges) 9.7352 14.2289 0.8671 0.9940 
Proposed method(Vertical and Horizontal edges) 4.8643 28.2702 0.7691 0.9931 
Proposed method(Common edges found between 13.41 5.563 0.9461 0.9957 
Horizontal and Vertical edges) 

Cup Image Canny 28.3816 0.44323 0.8133 0.9930 

(3264 x Sobel 28.44 3.0576 0.9648 0.9959 

1836) Proposed method(Vertical edges) 11.134 9.9401 0.9061 0.9939 
Proposed method(Horizontal edges) 11.4286 9.1913 0.9120 0.9940 
Proposed method(Vertical and Horizontal edges) 7.1055 17.5199 0.8496 0.9932 
Proposed method(Common edges found between 16.6041 2.6109 0.9687 0.9962 
Horizontal and Vertical edges) 

Parrot Image Canny 13.489 1.6127 0.8884 0.9957 

(1600 x Sobel 13.5014 4.6215 0.9559 0.9977 

1000) Proposed method(Vertical edges) 7.631 35.2328 0.7370 0.9939 
Proposed method(Horizontal edges) 8.1568 28.255 0.7759 0.9943 
Proposed method(Vertical and Horizontal edges) 3.4036 60.226 0.6222 0.9932 
Proposed method(Common edges found between — 10.5575 11.6803 0.8917 0.9960 


Horizontal and Vertical edges) 


4. CONCLUSION 

A novel approach for distinguish the edges for a given grey image is proposed. Edges are found 
along with noises for the lower resolution image. The results are improved over the existing approach for 
high resolution image. The proposed method uses the histogram technique to distinguish the edges within the 
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image. The algorithm utilizes the discontinuity of the intensity distribution to mark the edge information. 
Even the algorithm is tested with Gaussian noise image and salt & pepper noise image, the results are better 
when compared with the existing approach. At a certain noisy type of images, edges of the texture are more 
visible than by the other methods. The depiction of the depth is also visible in the output. The algorithm can 
be enhanced by adjusting the threshold value to detect the prominent edges. The edge can be traced with the 
width of one pixel on further enhancement. Since it can trace the edges for every intensity bin of the given 
images, this will be more suitable for medical image analysis. The major advantage of the proposed work is 
rejection of the false edges for the nearer coherent region which sometime perceived to be two different 
region. This may be helpful in analyzing the medical images and other edge detection application. 


REFERENCES 

[1] Ruigi Guo, Qieyun Dai and Derek Hoiem, “Single Image shadow Detection and Removal using Paired Regions”, 
International Conference on Computer Vision and Pattern Recognition, 2011, pp. 2033-2040 

[2] Lalit Gupta and Sukhendu Das, “Texture Edge Detection using Multi-resolution Features and SOM”, Pattern 
Recognition, 18th International Conference on Pattern Recognition, Vol. 2, 2006, pp. 199-202. 

[3] Ao Huanhuan and Yu Nenghai, “Edge Saliency Map Detection with Texture Suppression”, Image and Graphics 
(ICIG), 6th International Conference on Image and Graphics, 2011, pp. 309-313. 

[4] Sankar Krishnamurthy, S. Sitharama Iyengar, Ronald J. Holyer and Matthew Lybanon, “Histogram Based 
Morphological Edge Detector”, JEEE Transactions on Geoscience and Remote Sensing, Vol. 32, No.4, , July 1994, 
pp. 759-767. 

[5] Byung-Gyu Kim and Dong-Jo Park, “Novel Non contrast-Based Edge Descriptor for Image Segmentation”, JEEE 
Transactions on Circuits and Systems for Video Technology, Vol. 16, No.9, September 2006, pp. 1086-1095. 

[6] Qian Xu, Srenivas Varadarajan, Chaitali Chakrabarti and Lina J. Karam, “A Distributed Canny Edge Detector: 
Algorithm FPGA Implementation’, JEEE Transactions on Image Processing, Vol. 23, No. 7, July 2014, pp. 2944- 
2960. 

[7] R. Deriche, “Using canny criteria to derive a recursively implemented optimal edge detector’, International 
Journal of Computer Vision, Vol. 1, no. 2, 1987, pp. 167-187. 

[8] L. Torres, M. Robert, E. Bourennane, M. Paindavoine, “Implementation of a recursive real time edge detector using 
retiming technique”, Proc. Asia South Pacific IFIP International Conference on Very Large Scale Integr., 1995, 
pp. 811-816. 

[9] F. G. Lorca, L. Kessal, D. Demigny, “Efficient ASIC and FPGA implementation of IIR filters for real time edge 
detection”, Proc. IEEE ICIP, Vol. 2, October 1997, pp. 406-409. 

[10] J.F. Canny, “Finding edges and lines in images”, M.I.T. Artificial Intelligent Lab, Cambridge, MA, TR-720, 1983. 

[11] Sergio Carrato and Giovanni Ramponi, “Edge Detection Using Generalised Higher Order Statistics”, Higher Order 
Statistics, IEEE Signal Processing Workshop on Higher Order Statistics, 1993 proceedings, pp. 66-70. 

[12] Vishal Jain, Benjamin B. Kimia and Joseph L. Mundy, “Background Modelling Based on Subpixel Edges”, [EEE 
International Conference on Image Processing, Vol. 6, 2007, pp. VI-321 — VI-324. 

[13] Sangwook Lee, Sanghun Lee and Chulhee Lee, “Unsupervised Segmentation of Hyperspectral Images based on 
Dominant Edges”, International Conference on Computer Vision Theory and Applications (VISAPP), Vol. 1, 2014, 
pp. 588-592. 

[14] A.Leone and C.Distante, “Shadow detection of moving objects based on texture analysis’, Pattern Recognition, 
Vol. 40, Issue 4, April 2007, pp. 1222-1233. 

[15] Vidya.P, S.Veni and Narayanankutty K.A., “Performance Analysis of Edge Detection Methods on Hexagonal 
Sampling Grid’, International Journal of Electronic Engineering Research, Vol. 1, No. 4, 2009, pp. 313-328. 

[16] Abhishek S.N, Shriram K Vasudevan, Sundaram R.M.D, “An enhanced and efficient algorithm for faster, better 
and accurate edge detection”, Communications in Computer and Information Science, Vol. 679, 2016, pp. 24-41. 

[17] Sowmya V, N. Mohan, Soman K.P, “Edge detection using sparse banded filter matrices”, Second International 
Symposium on Computer Vision and the Internet (VisionNet’15). Elsevier Procedia Computer Science Journal, 
SCMS School of Engineering, Aluva, Kochi, 2015, pp. 10-17. 

[18] John Canny, “A Computational Approach to Edge Detection’, JEEE Transactions on Pattern Recognition and 
Machine Intelligence, Vol. PAMI-8, No. 6, Nov 1986, pp. 679-698. 

[19] ESA/Hubble & NASA, “Spiral Galaxy NGC 6384”, photographed by Hubble Space Telescope, February, 2011. 

[20] “Gorilla” image,http://gigglepedia.com/gigglepedia/p/animal-wallpapers.asp?id=1 &s=84 

[21] “Monkey” image, http://zibrato.com/cute-baby-monkey-background-wallpapers/ 

[22] Muammer Catak, “A NonLinear Directional Derivative Scheme for Edge Detection’, International Journal of 
Electrical and Computer Engineering (IJECE), Vol. 2, No. 4, August 2012, pp. 563-570. 

[23] ZhuZhong Yang, Jiliu Zhou, FangNian Lang, “Detection Algorithm of Airport Runway in Remote Sensing 
Images”, TELKOMNIKA Indonesian Journal of Electrical Engineering, Vol.12, No.4, April 2014, pp. 2776-2783. 

[24] R. Vijaya Kumar Reddy, K. Prudvi Raju, M. Jogendra Kumar, L. Ravi Kumar, P Ravi Prakash, S. Sai Kumar, 
“Comparative Analysis of Common Edge Detection Algorithms using Pre-processing Technique”, /nternational 
Journal of Electrical and Computer Engineering (IJECE), Vol. 7, No. 5, October 2017, pp. 2574-2580. 

[25] P. A. Khaire and Dr. N. V. Thakur, "A Fuzzy Set Approach for Edge Detection", International Journal of Image 
Processing (IJIP), Volume 6, Issue 6, 2012, pp. 403-412. 


Edge Detection Using Histogram Localization (Mohamed Asharudeen J) 


